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Abstract 

In this article we study the behavior of a group of economic agents in the context of cooperative 
game theory, interacting according to rules based on the Potts Model with suitable modifications. 
Each agent can be thought of as belonging to a chain, where agents can only interact with their 
nearest neighbors (periodic boundary conditions are imposed). Each agent can invest an amount 
(j£ = 0, ...,<? — 1. Using the transfer matrix method we study analytically, among other things, the 
behavior of the investment as a function of a control parameter (denoted (5) for the cases q = 2 
and 3. For q > 3 numerical evaluation of eigenvalues and high precision numerical derivatives are 
used in order to assess this information. 
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I. INTRODUCTION 



Consider a game where N players can invest their money (up to some upper limit) on 
some public-fund asset. The fund manager, as a rule, doubles the amount of money received 
and divides it equally among all N investors. Depending on the total amount each player 
invested, some might end up making a profit while others may lose money. Assume that 
players have no information whatsoever about their co-players' moves. The question is: 
what is the best move a player can make? If we adopt one of the tenets of classical game 
theory, namely that players are completely rational, then there are two possible solutions 
to the problem which maximize profit (when all invest the maximum amount possible, thus 
doubling their initial capital) or minimize losses (no one invests anything). In the real 
world however people are not rational in the sense of classical game theory and factors as 
expectations about the behavior of other players or some sort of insider information might 
play a role when deciding how to invest. 

The irrationality of market agents is one of a myriad of factors which account for the 
high complexity of financial markets and the difficulty in modelling them. Markets may also 
be affected by political turmoil, unseasonable weather variations and the like. In the past 
few years models have been introduced in order to throw some light into the behavior of 
markets, mainly with aims at forecasting long-term behavior see for example [3,0]. These 
models, which have the advantage of being either analytically or numerically treatable and 
usually use some kind of data input from real markets, are nonetheless unable to take into 
account the human factor in decision-making scenarios. To circumvent these difficulties a 
new approach, inspired on the ideias of statistical mechanics has been suggested 0], where 
one extends the set of causal factors in decision-making scenarios from the individual- 
specific to group determinants of behavior: players' decisions are not market-mediated but 
rely on group-level influences. 

With these ideias in mind our aim in this work is to extend the model for the game 
discussed above through the introduction of cooperation, i. e. we allow agents to have partial 
information on the decision of its immediate neighbors, in a way to be described precisely in 
what follows. Furthermore, we allow for some kind of randomness, the only thing known a 
priori being the probability of some decision, and not the decision itself. In this way we hope 
to describe the average behavior of a large group of agents without entering in the details 
of a realistic (and certainly very difficult) theory on psychological state of each agent. 

Our main interest will be to see how is the average behavior of a (large) group of cooper- 
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ative economic agents. Each agent, labelled by the index i, is allowed to invest an amount 
<jj. Before the investment, the agent i exchange information with agent i + 1 (defined as the 
neighbor of i; this concept is symmetric, i.e., i is also neighbor of i + 1). Based on that infor- 
mation agents make decisions as to how much they will invest according to some probability 
distribution parameterized in terms of a two real variables: J and (3. The first measures how 
strongly people interact with each other (group-level influence) while the later is a measure 
of how strongly a player might deviate from the group. To model this we choose a suitably 
defined function which measures the probability of i investing <jj given that i + 1 would like 
to invest (J i+ i. In a way to be precisely formalized later, j3 allows us to change the expected 
behavior of each agent. 

The paper is organized as follows: In section 2 we explain the model and make the 
connection with statistical mechanics. Using the standard transfer matrix technique [g] we 
analyze in section 3 some integrable cases (q = 2 and q = 3) in order to gain information 
about how the average investment changes as a function of (3. In section 4 we numerically 
evaluate the evolution for q > 4. We introduce a method for calculating the investment using 
derivatives of the biggest eigenvalue, which is based on the use of 5 points in a graphic. We 
finish the paper with some conclusions and perspectives. 

II. FORMULATING THE PROBLEM: THE POTTS MODEL 

We consider an ensemble of N agents, where each can invest an amount oi G {0, 1, . . . , q— 
1}, q a fixed integer. This restriction on q has been made for the sake of clarity. The 
methods employed can be easily generalized to the case where <7j G {do, . . . , d q _i}, the c^'s 
being arbitrary real numbers. 

The families of conditional probabilities, i. e. the probability that i + 1 would invest a i+ i 
given that i invested cij are chosen as 



Our motivation for this particular choice comes from physics, where this probability is 
interpreted as that of two variables (called classical spins) <jj and o"j + i being equal or having 
different values. It depends on two physical parameters: J(o"j) is the so-called interaction 
strength. This is in most cases a material-dependent quantity and accounts for the different 
collective properties materials may exhibit (ferromagnetic, antiferromagnetic, etc.). (3 is 
proportional to the inverse temperature T -1 and brings about entropic effects (ordered 
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states for low temperatures and disorder for high temperatures). Being proportional to 
a temperature, in physical systems (3 is always non negative, and we likewise assume our 
> 0. 

In general J and f3 can be seen as competing terms: J is associated to the energy cost of 
a given spin configuration. For J > (< 0) a configuration where spins are equal (different) 
has a higher energy than the opposite configuration, which means that it is energetically 
more favorable to be non-magnetic (or ferromagnetic); on the other hand, the temperature 
P^ 1 tends to destroy magnetic order. Transposing these ideas to the financial context can 
say that J((Xj) measures how strongly people respond to their neighbors' moves, that is if 
they are susceptible to the influence of other players or not. In this sense it models distinct 
profiles of investors, which can go from agressive (does not go along "with the pack") to 
conservative (does what others do). (3 is a measure of the strength of individual response and 
independent of what others do. Social scientists refer to this term as the "individual-specific 
random" and "unobservable" (from the point of view of the modeler) since it is associated 
to personal beliefs ^]. 

With the neighbor-to-neighbor interaction rule introduced above we can describe the 
behavior of the whole group: The first important quantity which needs to be defined is the 
joint probability density (j.p.d), P(o~i, <7 2 , Cat) of a particular investment configuration 
a = (cri, 02, crjv) of a group. The quantity invested is defined through 

L(a) = oi + a 2 + ••• + &N (2) 

and we would like to obtain the average value of L. 

To calculate the j.p.d we may adopt a recursive formulation without any loss of generality: 
we take the investment of the first agent to be exactly <Ti, i.e., P(o~i) = 1 and from that 
derive the quantity we want. With this "boundary condition" on P{ci) we have the following 
theorem: 

Theorem 1 The probability distribution P(o~i, o~ 2 , o~n) can be written as a product form 

N-l 

P(a 1 ,a 2 ,...,cr N )=Y[P( 'i+iWi) ( 3 ) 

i=i 

with P{<Ji\ai- U cri) = P(a i \a i - 1 ) for i = 1, N. 

Proof. /From the definition 

P(a N \a u cr 2 , 0at_i) = — (4) 

-r(01, 02, •••,0JV-lJ 



Considering the hypothesis P((7jv|ci, a 2 , ...,crjv-i) = P{pn\on-\) we thus have 



P(a 1 , a 2 , a N ) = P(a 1 , a 2 , o-Ar_i)P(crAr + i|a-Ar) . 



(5) 



and applying this recursively: 



P(a 1 ,a 2 ,...,cr N ) = P(a l ,a 2 , ...,a N _i)P(a N \a N _ l ) 



P(ai,a 2 , ...,a N - 2 )P(a N -i\a N - 2 )P(a N \a N - 



N-l N-l 

p(o-i) n ^(^i+iko = n 



According to equations ((T} and (jHJ) we have 



P(a 1: a 2 , ...,a N ) = c exp 



AT-1 



i=l 



where c is the normalization constant defined before and such that 

dq-l d q -l dq-l 

O"! =(2q CT2 =rfo crN=do 



A. Investment formulas and the Potts model 



The Potts hamiltonian of N interacting spins under the action of a magnetic field D is 
given by p 

N 

H(a) = J2A°i)^+Dj2 a *- ( 6 ) 

The fact that the total investment L (J2j) is mathematically the same as the magnetization 
of the Potts hamiltonian (jUJ) means that we may directly transpose the techniques and ideas 
of statistical mechanics into the financial scenario. 

The probability density function for the system to present a specific investment value L 
is given by 

' Zjexpl-PHi*)] if L = J2f =l a l 

P(L) = \ (7) 
otherwise 
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where Zn(/3), the normalization constant, is a sum over all possible configurations 

Zn(P)= ^Pi-PH(a)}. (8) 

(7j=0,...,q— 1 

and is known as the partition function. 

Our aim is to describe how the investment depends on /3, given a fixed set of parameters 
J(l), J(2), J(N). For this purpose, let us consider the expected value of L according to 
the distribution (JJJ), i.e. 

(L( ( x)) = l fe'ijexp 

o-j=0,...,g-l \i=l / 

For the sake of those not familiar with the methods of statistical mechanics, we briefly 
discuss how in our analogy between spin systems and economic games quantities of interest 
can be calculated: 

N 

1. The term D^ai is introduced for convenience since investment per capita is 

i=i 

calculated through the formulae 

l ^ = ^ L ^=-jNW l °^\ D __ a ' 

which is the analogue in statistical mechanics to the average spontaneous magnetiza- 
tion. It clearly obeys the inequality < I < q — 1. 

2. As previously discussed, investment depends essentially on the way how neighbors 
cooperate, i.e. the distribution J(0), J(l), .... J(q — 1) defines the kind of profiles 
investors have. 

The first question one might ask would be: what is the behavior of the per capita in- 
vestment l(f3) at /3 = and /3 — > oo ? One may describe these limits in a straightforward 
manner: 

Theorem 2 The per capita investment is such that 

/(oo) = a m in if Jmin satisfies the inequality J min < for all k = 0, . . . , min — 1, min + 

1,. 
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Proof. The probability of a given state a = {(Ji}f =l is 

For /3 — all states are equiprobable since, from the equation above, the probability does 
not depend on a; hence, the probability of each state is l/q N ■ For simplicity we assume that 
q N is even, the odd case being left to the reader. The per capita investment can be written 

as 



1(0) = — N (number of states with sum = i)i 



N(q-l) 
i=0 

We can see that the number of states corresponding to the sum i is the same as the number 
corresponding to the sum N(q — 1) — i. Then, 



AT(g-l)/2 

1(0) = N (number of states with sum = i)N(q — 1) 



Nq 1 

y i=0 



A^g w y 2 2 

since the number of states with sum between and N(q — l)/2 correspond to (q N )/2. 
For arbitrary values of (3 let us call cr M = {a- min }^ =1 . Then, the probability of any state 

is 

P(a) 



e -f3 H(a') 



In the limit (5 — > oo, -P(c) = except for the state <7m, where P((Jm) = 1- ■ Notice that 
the case where minj Jj is reached in more than one point is not covered by the theorem. 

In the next section we give a case by case description of investment as a function of f3. 
To do this we adapt transfer matrix method to our model. Contrary to spin systems, in 
the present problem nontrivial behavior appears, and this might lead to interesting new 
possibilities in the scenario of economic games. 
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III. ANALYTICAL CASES 



A. The two state model: q = 2 (Ising Model) 



We start by considering the partition function 



(71=0,1 a N =0,l 



N 



N 



Z = E ■■• E e x P[-PJ2J(< T i) S wn.i-2 D 'E( <T i +(T i+l. 



i=l 



i=l 



A? 



E ••• E nexp(-^(a j )^ j+1 -f% + <7 W )) 

0-1=0,1 <7JV=0,ll=l 



(9) 



M 



= Tr M N 

where <Jn+i = &\ (periodic boundary conditions are assumed) and the transfer matrix M is 
given by 

exp(-/?J ) exp(-((3/2)D) 
exp (-((3/2)D) exp (-pj x - (3D) 

It is possible to diagonalize M and write Q in terms of eigenvalues Ai and A 2 of M in a 
simple way 

Z N (P,D) = \» + A 2 V (10) 

These eigenvalues are given by 



with 



A li2 (/5) = i[(e-^+e- D ^)±v / e" 



_ 4 e -^/3 _|_ e -2G9Jo) + e -2/9(D+Ji) _ 2e -P(D+Jo+Ji) 

J, = J(i) i = 0,l . 



One can now write 



log(Af+A?)| 



i a 

/9AT ao v^i 1 "2 ) |D=0 



i i 

/3V (X^+W) 



D=0 



11 1 [ d \ | (\ 2 \N-1 d \ l 

5A a (l+CAa/Ai)™) ISD Al Ui^ aD A2 J 



D=0 
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In the limit of a large group of agents (N — > oo) the expression above converges to 



D=0 



An explicit calculation gives 
l(/3, J , Ji) 

where 



_ e /3(Ji-Jo) + 2e 2 ^ + 1 + y/A 



4 e 2/3Ji _ 2 e P(Ji-Jo) + e 2/3(Ji-j ) _|_ e P{Ji-Jo)y/A + ^/A + 1 



A(/3, Jo, Ji) = 4e 2/3Jl - 2e /3Jl e-^° + e - 2/3Jo e^ Jl + 1 . 
A few remarks can be drawn from these equations: 

1. With the explicit expressions of the eigenvalues one can easily show that 

1{P = o) = I 

independently of the values of Jo and J\ (as expected from the last theorem). 

2. When J = J\ one has A = 4e 2/3Jo and 



(11) 



(12) 



l{(3) = 1/2 for all (3 E R 

The other important limit l(/3 — > oo) has an explicit Jj dependence that can be summarized 
below: 

\ if J > Ji > or Ji > J > 

l((3 oo) = < 1 if J x < J < or Ji < and J > 

k if J < Ji < or Ji > and J < 
This result follows from the expression for l((3, J , J\) above. At this point we would 
like to make an important comment on the mean value of Jj. More specifically we are 
interested in the positivity or negativity of Jj for this tells us to what extent agents are 
cooperating and the probability of their making investments. For example, if Jo > and 
Ji < 0, the probability of two neighbors investing the same <jj = 1 is greater than the 
co joint investment of <j{ = 0. 

This behavior can be better seen in Figs. Q and El The first one depicts investment as a 
function of (3 for the ratio Jq/J\ > 0. In the second figure the sign of this ratio is reversed. 



9 



> 




FIG. 1: Investment as function of (3. Case Jq/Ji > 0, q = 2. 

These clearly show how investment behavior (how agents cooperate) is drastically modified 
as a function of the profit and depends not only on cash flow (cash receipts minus cash 
payments over a given period of time). 



50.4 




FIG. 2: Investment as function of (3 to q = 2, case Jq/Ji < 
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B. The three state model: q = 3 



In this case we have the following transfer matrix: 



M 



,-/3D/2 



,-BD 



\ 



e -3D/2 e -/3(J 1+ D) e -3f3D/2 
e -8D e -3/3D/2 e -8(J 2 +2D) 



The problem of computing eigenvalues is analytically tractable for general values of set 
{Jo, J±, J2} but the expressions obtained are generally very difficult and do not improve our 
understanding of the problem. However some interesting sub-cases can be considered: 

1. Case 1: Jq = J\ = and J2 = J 

In this case by making the change of variables 

x = e~P J and y = e~^ 2 



we arrive at 



M 



1 ' 1 y y 2 ^ 

y y 2 y 3 
\y 2 y 3 xy A ) 



(13) 



One may clearly see that in this matrix the second row is obtained from the first through 
multiplication by y. Hence detM = and therefore M admits as an eigenvalue. One may 
thus compute analytically the other eigenvalues by solving the equation 

A 2 - (1 + y 2 + xy 4 )X + (x - l)(y 4 + y 6 ) = 0, 



where the largest eigenvalue is 



A = ^(l + 2 / 2 + a;|/ 4 + v / A) 



Here we have 



A = 2y + 5y 4 + Ay b - 2xy A - 2xy b + xY + 1 =>- A(D = 0) = 12 - Ax + x 
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and so 



1 8\ dy 
f3X dy dD 



D=0 



(2+x+y/12-4x+x 2 ) 



1 + 2x 



12-5x+2x 2 
Vl2~4x+x 2 



I _|_ 2e~ l3J + 12 - 5e ~ /3J + 2e ~ 2 ' 3 ' 7 

v /l2-4e-/ 3J +e- 2 / 3 -' 



( 2 +e- l '3^ + v /l2-4e-/ 3 J+e- 2 / 3J ) 

One may observe that two cases follow from theorem 3.1: (3 = and /3 — > oo when J < 0. 
We have Z(0) = 1 and if J < 0, i(/3 — > oo) = 2. But with the expression for / we can also 
obtain results beyond the range of the theorem, for example, when J > 0. In this case 
l{(3 -»• oo) = [1 + Vl2] = 0.816 97- ••. 



£ Case £: J = J2 = and J a = J 



In this case we have 



M 



/ 1 y y 2 ^ 

y xy 2 y 3 



2 3 4 

y y y 



The discussion is analogous to that of section (fill B 1|) . The largest eingenvalue is given by 

1 



A = - 
2 



(y 4 + xy 2 + 1) + VA 



with 



A = Ay 2 + 2y 4 + % 6 + y 8 - 2xy 2 - 2xy 6 + x 2 y A + 1 . 
A straighforward calculation gives 

l(/3) = 1 

Thus the investment is independent of the parameter J x when J = J 2 = 0. 



3. Case 3: J = J and J\ = J 2 = 



Now M takes the form 



M 



^ x y y 2 * 

y y 2 y 3 

2 a 4 
y y y 



12 



As before we have 



a = I 

2 



x + y 2 + y 4 + VA 



with 



A = x 2 + Ay 2 + 5y 4 + 2y 6 + y 8 - 2xy 2 - 2xy 4 = x 2 + 12 - 2x - 2x = 12 + x 2 



Ax 



For we have the following expression 

1 



w 



3 + 



(12 - 3e-^ J ) 



V12 - Ae-P J + e" 2 ^ 



e -(3J + 2 + v / 12-4e-^ J + e- 2 ^ J 
The case J > is not covered by theorem 3.1. From the expression above we obtain 

... . 3 + ^12 

The only cases where one has analytical solutions are for q < A. For other values of q one 
has to employ numerical methods in order to gain some information, as we discuss in the 
next section. 



IV. NUMERICAL ANALYSIS 



For those cases which are not analytically treatable we can employ an algorithm that 
combines a routine of numerical derivation with eigenvalues computing. To see how the 
method work, we first consider the following matrix, written as a function of £: 



/ 



M(0 



-f3 Jo 



-(3 €/2 



e -^/2 e -(3 (Ji+O e -3/3£/2 

e -K e -mm e -/3 {J2+20 



-/3£(«-l)/2 e -fcq/2 e -p£(q+l)/2 



e -tH{q-l)/2) \ 
e -^(g+l)/2 



-(3 (J,_i+(?-l)0 



Let Ag be the largest eigenvalue corresponding to matrix M(£) and A_,e the one from M(— £). 
If £ << 1 one has: 



D=0 



Ag - A_^ ~ e (2fc+1) A (2fc+1) (0) 
2£ 2(2* + 1)! 



since 



A, = Ao + A' £ + ^ AS + Ja " + 



(14) 



A — A - A' £ + ^-Aq - li-A'o" + .. 
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A numerical estimate to order ( 0(£ 2 ) ) is 

A"»A({) = 

A more refined numerical estimate can be obtained using not only two but four points Ag, 
A-£, A2^ and A_2§. 
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FIG. 3: Investment as function of (5 to q = 3. A comparison of the numerical result with analitical 
results. 



Theorem 3 Consider a function A G C c 
of \'(D) is given by 



A numerical approximation to order 0(£ 4 



A( 2 )A(0 = -iA«A(20 + ^A( 1 )(0 



Proof. Considering a Taylor expansion 



A 2? = A + 2£A' + 4£ 2 | + 8f + 16^ + 



4! 



(15) 



and 



i An „ A 



(3) 



A 



(4) 



A 2€ = A - 2£A' + 4£ 2 ^ - Sf 5 -^ + 16^-^ + •••• 



Combining the equations (JHJ), ([To]) and (fTfjj) . we obtain 



(16) 



5\(5) 



- A 25 + A_ 2€ - 8A_ 5 = 12£A' - -e 5 A 



So 



A' 



A 2 ^ — A_ 2 f 


4 




L ^ 


~ 3 


2£ 



14 



which gives us 

A (2) A(0 = -V)A(2£) + ^A«(0 

■ 

To assess the applicability and performance of the method, we applied it to the integrable 
q = 3 case with numerical results using 4-point derivative. In Fig. 3 we show how the 
numerical results compare with the analytical ones. By using a higher number of points we 
observe a significant difference on the results (see Fig. 4) over selected regions, as compared 
to a lesser number. 

A. Numerical analysis for q > 3 considering distinct profiles of agents 

In this section we analyze some numerical results for q > 3. We consider three possible 
profiles: 

1. aggressive or risk-prone agents: in this case the probability of an agent's invest- 
ment increase as a function of invested value. We modelled this through 

J(a t ) = -(a t + 1) < 0, 

where <7j — 0,...,q — 1. 

2. conservative agents: the probability of investment decreases as function of invested 
value. In this situation 

where <7j — 0,...,q — 1. 

3. random agents: the probability of the investment is randomly chosen for each agent, 
such that 

J{a t ) = Lrand[0, 1] • q\ 
where rand[0, 1] is a random number uniformly generated in the interval [0,1]. 

We generated plots with q = 10, 15, 20 for three different profiles. In Fig. 5, we show 
the risk-prone (a) and conservative (b) profiles. From (a) we conclude that all agents are 
inclined to invest the maximum quantity for j3 — * oo since greater quantities are privileged 
by the probability distribution. Differently, for (b) as (3 — > oo , the investment of each agent 
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FIG. 4: Comparative numerical derivates for q = 3, using three points and five points, (a) Jo = 
J\ = and J2 = J; (b) Ji = J2 = and Jo = J; (c) Jo = J2 = and J\ = J. 

goes to 0. We then may conclude that conservative agents lead to the situation of complete 
stagnation as (3 — ► 00 , independent of the number possibilities in the investment q. On the 
other hand, risk prone agents lead the market to invest the maximum at this limit. 

An alternative profile seems to be more appropriate: The random profile (3) was also 
explored in an experiment using 12 seeds (12 different random choices of the string J(cjj), 
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FIG. 5: Numerical results for the behavior of investiment as function of (3 for q = 10, q = 15 and 
g = 20 for two diferent profiles: (a) agressive and (b) conservative. 

% = 1...N) and q = 15 as depicted in Fig. 6 ). This figure represents the average over the 
12 seeds. The behavior of seeds are not similar in the sense that they may yield different 
values of investment at (3 — * oo. 



V. SUMMARY AND CONCLUSIONS 



We have studied the investment behavior of a group of agents as a function of a parameter 
that mimics the mean profit obtained by agents. Our results illustrated different situations 
based on possible investors' profiles. In a model where q represents the number of possible 
investment amounts, we obtained analytical results for q = 2 and 3. For larger values of q we 
performed a series of numerical simulations by combining exact diagonalization algorithms 
with numerical derivatives. 

Our results indicate that the behavior of each investor is key to determining the dy- 
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FIG. 6: (a) Plot for 12 different seeds for q = 15; (b) mean value over the seeds. 

namics of the market. As recent results in the context of agents' simulation show [1], pure 
mathematical models can capture some of the intrincacies of real markets when they, as 
pointed out in try to include real people's idiosyncrasies (beliefs, sentiments, etc.) that 
are known to play a significant role (not to mention other important influences as season- 
able changes in production, political turmoil, and the like). Even though our model is still 
mathematical, in the sense that it is based on a well known model of statistical mechanics 
and we identify behavior in terms of known physical quantities, we believe that our ideas 
might help indicate a way towards a more realistic market modelling. 
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